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ABSTRACT 

We present Westerbork Synthesis Radio Telescope Hi images, Lovell telescope 
multibeam Hi wide-field mapping, William Hcrschcl Telescope longslit cchcllc Can 
observations, Wisconsin Ha Mapper (WHAM) facility images, and IRAS ISSA 60 
and 100 micron coadded images towards the intermediate velocity cloud (IVC) at 
' +70 kms -1 , located in the general direction of the M 15 globular cluster. When corn- 

er . bined with previously-published Arecibo data, the Hi gas in the IVC is found to be 

CO ■ clumpy, with a peak H I column density of ~ 1.5 x 10 20 cm -2 , inferred volume density 

(assuming spherical symmetry) of ~ 24 cm -3 / D (kpc), and maximum brightness 
temperature at a resolution of 81" x 14" of ~ 14 K. The major axis of this part of the 
IVC lies approximately parallel with the Galactic plane, as does the low velocity Hi 
gas and IRAS emission. The Hi gas in the cloud is warm, with a minimum value of the 
• full width half maximum (FWHM) velocity width of 5 kms" 1 corresponding to a ki- 

netic temperature, in the absence of turbulence, of ~ 540 K. From the Hi data, there 
are indications of two-component velocity structure. Similarly, the Can spectra, of 
' resolution 7 kms -1 , also show tentative evidence of velocity structure, perhaps indica- 

^ \ tive of cloudlets. Assuming there are no unresolved narrow-velocity components, the 

mean values of logi (N(Caii K) cm -2 ) ~ 12.0 and Call/Hi ~ 2.5 xlO -8 are typical of 
observations of high Galactic latitude clouds. This compares with a value of Can/Hi 
> 10 -6 for IVC absorption towards HD 203664, a halo star of distance 3 kpc, some 
3.1 degrees from the main M 15 IVC condensation. The main IVC condensation is 
' detected by WHAM in Ha with central LSR velocities of ~ 60—70 kms -1 , and inten- 

sities uncorrected for Galactic extinction of up to 1.3 Rayleigh, indicating that the gas 
is partially ionised. The FWHM values of the Ha IVC component, at a resolution of 
1 degree, exceed 30 kms -1 . This is some 10 kms -1 larger than the corresponding Hi 
value at similar resolution, and indicates that the two components may not be mixed. 
However, the spatial and velocity coincidence of the Ha and Hi peaks in emission 
towards the main IVC component is qualitatively good. If the Ha emission is caused 
solely by photoionisation, the Lyman continuum flux towards the main IVC condensa- 
tion is ~ 2.7xl0 6 photons cm -2 s -1 . There is not a corresponding IVC Ha detection 
towards the halo star HD 203664 at velocities exceeding ~ 60 kms -1 . Finally, both 
the 60 and 100 micron IRAS images show spatial coincidence, over a 0.675° xO. 625° 
field, with both low and intermediate velocity H I gas (previously observed with the 
Arecibo telescope), indicating that the IVC may contain dust. Both the Ha and ten- 
tative IRAS detections discriminate this IVC from High Velocity Clouds although the 
Hi properties do not. When combined with the Hi and optical results, these data 
point to a Galactic origin for at least parts of this IVC. 

Key words: ISM: general - ISM: clouds - ISM: individual objects: Complex gp - 
ISM: structure - globular clusters: individual: M 15 - radio lines: ISM. 
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1 INTRODUCTION 

The study of intermediate velocity clouds (IVCs) remains 
one of the most challenging in contemporary Galactic as- 
tronomy, with several issues concerning IVCs remaining un- 
resolved. These include, but are not limited to, the method 
of their formation, their relationship (if any) with high veloc- 
ity clouds (HVCs), and the question as to whether IVCs are 
sites of star formation in the halo of the Galaxy (Kuntz & 
Danly 1996; Christodoulou, Tohline & Keenan 1997; Ivezic 
& Christodoulou 1997). This latter possibility is under- 
pinned by the fact that within the Galactic halo, there exists 
a population of early B-type stars whose velocities, ages and 
distances from the Galactic plane (z) are incompatible with 
them being formed within the disc. A possible site for their 
formation is IVCs/HVCs via cloud-cloud collisions and sub- 
sequent compression of the gas (Dyson & Hartquist 1983). 
Such collisions are thought to be a viable star formation 
mechanism within at least the discs of galaxies, albeit where 
the gas density and cloud-cloud collision rates are somewhat 
higher than inferred in IVCs/HVCs (Tan 2000). 

The solution to both the star formation question and 
also any possible relationship between HVCs and IVCs re- 
quires both the analysis of aggregate parameters of well- 
defined samples of IVCs and HVCs, and also more detailed 
studies of individual objects. In this paper we report on radio 
Hi aperture synthesis, Hi multibeam wide field mapping, 
longslit Call observations, Wisconsin Ha Mapper (WHAM) 
facility images, and IRAS sky-survey archive data retrieval 
towards a particular IVC located in the general direction of 
the M15 globular cluster (RA=21 h 29 m 58.29 s , Dec=+12° 
10' 00.5" (J2000); Z=65.01°, 6=-27.31°). These observa- 
tions are amongst the first Hi synthesis data to be taken 
of positive-velocity IVCs, which remain poorly-studied as a 
group of objects. 

The M 15 H I cloud lies at a velocity of ~ +70 kms -1 in 
the dynamical Local Standard of Rest (Allen 1973) ; its dis- 
tance tentatively lies between ~ 0.8-3 kpc (Little et al. 1994; 
Smoker et al. 2001a). The upper distance limit is gleaned 
from the fact that IVC absorption at ~+70 to +80 kms -1 
is observed in the spectrum of HD 203664, a halo star of 
distance 3 kpc and ~ 3.1 degrees from M 15 (Little et al. 
1994), combined with the detection of IV Hi approximately 
mid-way between the M 15 IVC and HD 203664. 

The deviation velocity of the M 15 IVC at such a mid- 
Galactic latitude puts it in on the borderline between the 
normal definitions for intermediate and high velocity clouds 
(c.f. Fig. 1 of Wakker 1991), although in common with Sem- 
bach (1995) and Kennedy et al. (1998) here we classify it 
as an IVC. The line-of-sight position of the M 15 IVC is be- 
tween the negative-velocity Local Group barycentre cloud 
Complex G and the Galactic centre clouds (Fig. 8 of Blitz 
et al. 1999), hence the M 15 IVC is a part of IVC Complex 
gp (Wakker 2001). 

Previous observations in Hi emission using the Lovell 
and Arecibo telescopes (Kennedy et al. 1998; Smoker et al. 
2001a) have shown that the IVC consists of several conden- 
sations of gas spread out over an area of more than 3 square 
degrees, with structure existing down to the previous res- 
olution limit of ~ 3 arcmin. The brightest component is 
located towards M 15 itself and has peak H I column density 
at the Arecibo resolution of ~ 8xl0 19 cm -2 . In this paper, 



we study this part of the IVC at higher spatial resolution. 
The mass of this particular clump is ~ 20 D 2 Mq, (where 
D is the distance in kpc), thus for this particular object, in 
the absence of H2, there is insufficient neutral gas to form an 
early-type star. Low-resolution absorption-line Call and Nal 
spectroscopy (Lehner et al. 1999) towards cluster stars ten- 
tatively found cloud structure (or variations in the relative 
abundance) over scales as small as a few arcsec, with fibre- 
optic array mapping in the Na I D absorption lines (Meyer & 
Lauroesch 1999) obtaining similar results with structure vis- 
ible on scales of ~ 4 arcsec (velocity resolution ~ 16 kms -1 ). 
Using empirical relationships between the sodium and hy- 
drogen column densities, Meyer & Lauroesch (1999) derived 
values of the Hi column density towards the cluster centre 
of « 5xl0 20 cm -2 , some 6 times higher than that found 
using the Arecibo Hi data alone; the difference may be at- 
tributable to fine-scale cloud structure. Assuming spherical 
symmetry, a volume density of « 1000 cm -3 is implied by 
these latter results, similar to values obtained for gas in the 
local ISM (e.g. Faison et al. 1998, although see Lauroesch, 
Meyer & Blades 2000). Such a high volume density and im- 
plied overpressure with respect to the ISM perhaps indicates 
that the assumption of spherical symmetry is invalid and 
that there may be some sheet like geometry in the IVC as 
has been postulated for low- velocity gas (Heiles 1997). 

In the current paper, we extend our studies of the M 15 
IVC to higher resolution and different wavelength regions in 
order to investigate three areas. Firstly, Hi synthesis map- 
ping, WHAM Ha and IRAS ISSA survey data retrieval to- 
wards the IVC were performed to see if the H I, Ho and 
infrared properties of the M 15 IVC are compatible with ei- 
ther low velocity gas or HVCs in general, and whether there 
are any differences between the types of object, perhaps at- 
tributable to differences in the formation mechanism or cur- 
rent environment (for example, distance from the ionising 
field of the Galaxy). Secondly, wide- field medium-resolution 
H 1 data were obtained in search of more IVC components, to 
trace how the gas kinetic temperature changes with sky po- 
sition, and possibly determine the relative distance of cloud 
components from the Galactic plane (c.f. Lehner 2000). 
Thirdly, longslit echelle Call observations were undertaken, 
using the centre of M 15 as a background continuum source, 
in order to look for small-scale velocity and column density 
substructure within the IVC which could indicate the pres- 
ence of cloudlets, collisions between which in certain IVCs 
may be responsible for star formation in the Galactic halo. 

Section^ describes the observations and data reduction, 
Sect. gives the results, Sect. ^ contains the discussion and 
Sect. H presents a summary and the conclusions. 

2 OBSERVATIONS AND DATA REDUCTION 

2.1 Westerbork aperture synthesis Hi 
observations 

21-cm aperture synthesis Hi observations of the M 15 IVC 
were obtained during two observing sessions, each of 12 
hours, using the Westerbork Synthesis Radio Telescope 
(WSRT). The first in 19 December 1998 had a minimum 
antenna spacing of 32-m, the second in 17 April 1999 had 
a corresponding separation of 72-m. The velocity resolution 
for all observations was 1.03 kms -1 . 
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Standard methods within AlPSy were used to reduce the 
visibility data. Reduction included amplitude calibration us- 
ing 3C 286 and 3C 48 (assuming flux densities of 14.76 Jy 
and 15.98 Jy, respectively), phase calibration and flagging 
of bad data, and concatenation of the two UV datasets us- 
ing DBCON. The calibrated dataset was mapped with IMAGR 
using both quasi-natural and uniform weightings, with the 
Arecibo map from Smoker et al. (2001a) being used to set 
the locations of the CLEAN boxes for each velocity channel 
interactively. The respective beamsizes of the final images 
were 111" x 56" and 8l"xl4" (approximately North-South 
by East-West), with the corresponding rms noises being 2.3 
mjy beam -1 and 1.3 mjy beam -1 for the naturally and 
uniformly-weighted data. 

As the WSRT maps suffer from missing short-spacing 
information, it was decided to create 'total power' chan- 
nel maps by combining the current observations with the 
single dish Arecibo data. For this we used miriad (Sault, 
Teuben & Wright 1995). The two datasets were first regrid- 
ded so that they had the same central coordinate and chan- 
nel width. Following this, we multiplied the Arecibo data by 
the WSRT single-dish beam, converted to flux density units, 
and combined the WSRT and weighted Arecibo maps using 
immerge. Combination was performed by specifying an an- 
nulus of between 50-150-m in the Fourier domain where the 
high and low resolution images were made to agree. Within 
this annulus it was found that, where the Hi was strong, 
the flux density scales of the two datasets were the same. 
Hence the 'calibration factor' used when combining the two 
datasets was set to unity. The resulting combined channel 
maps were divided by the WSRT single dish beam to pro- 
duce the final combined images. These are of the resolution 
of the WSRT data (lll"x56"), but with the correct zero- 
spacing flux. Finally, moments analysis and Gaussian fitting 
were performed on the WSRT channel maps and combined 
cube from +50 to +90 kms -1 in order to determine the 
Hi column densities, velocities and velocity widths of the 
IVC. H I column densities were derived using the relationship 
Njj/=1.823xl0 18 / T B dv (cm -2 ), where v is the velocity in 
kms -1 . 



2.2 Lovell Telescope multibeam Hi observations 

Lovell telescope multibeam 21-cm Hi observations cover- 
ing a ~ 4.9x9.3 degree RA— Dec region, centred upon 
RA=21 h 22 m 47 s , Dec=+10°00'29" (J2000) were taken dur- 
ing 2 June 2000. The resolution of these data is 12 arcmin 
spatially and 3.0 kms -1 in velocity. The total integration 
time was 4 hours and the RMS noise per channel was ~ 
0.1 K. The observed field is not quite fully sampled, with 
the spacing between the beams being 10 arcmin, compared 
with the spatial resolution of 12 arcmin. During observing, 
data were calibrated in terms of antenna temperature by a 
frequently-fired noise diode. Off-line, bootstrapping of the 
flux scale to some previous Lovell telescope H I observations 
taken during 1996 enabled calibration in terms of bright- 
ness temperature (Tb). Data reduction procedures will be 



* AIPS is distributed by the National Radio Astronomy Observa- 
tory, U.S.A. 




Figure 1. Lovell telescope multibeam Hi spectrum towards 
RA=21 h 26 m 47 s , Dec=+09°16'23" (J2000). The dotted line 
shows the results of a single-component Gaussian fit. The bright- 
ness temperature calibration is accurate to ~ 40 per cent with 
the negative features present near to ~ kms -1 being an arte- 
fact of the baseline removal procedure. The dashed line is the 
(data— single Gaussian fit) residual. 

discussed in more detail in a future paper. For now, we sim- 
ply note that most of the process is automated, the bandpass 
removal being performed within AIPS ++ using livedata and 
the data being gridded using GRIDZILLA (Barnes 2001). As 
the bandpass is simply a running average of the observed 
spectra in a window, incorrect removal is performed around 
the low- velocity gas where there is no 'empty' region of the 
sky. However, at the velocity of the IVC (v ~ +70 kms -1 ) 
the baseline subtraction is adequate. An example of a spec- 
trum after automated bandpass removal is shown in Fig. 

After gridding, the data were imported into AIPS whence 
a column density map was made of the entire field bet ween 
+50 and +90 kms -1 which is presented in section 3.2. 



2.3 William Herschel Telescope longslit UES 
observations 

WHT Utrecht Echelle Spectrograph (UES) observations to- 
wards the centre of the M 15 globular cluster were taken 
during bright time on 19-21 July 2000. Twin EEV CCDs 
were used. In combination with the E79 grating, this gave 
near-complete wavelength coverage from 3800-6300 A, with 
a slitwidth of 0.85" providing an instrumental FWHM res- 
olution of ~ 6 kms -1 . On the first night of the run, a slit 
height of 16 arcseconds was used; this has the advantage 
that the full wavelength coverage is obtained with no over- 
lapping orders on either of the CCDs. On the second and 
third night's observations, in order to maximise the amount 
of sky available for the calcium lines at 3933. 66A, a slit of 
height 22 arcseconds was used which leads to overlapping 
orders in the red CCD. ThAr arcs were taken two or three 
times a night to act as the wavelength calibration with Tung- 
sten and sky flats also being obtained. The echelle proved 
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to be remarkably stable with the centres of the arc lines 
varying by less than 1.5 pixels throughout the whole of the 
run. 

The original intention of the observations was to ob- 
tain 7 parallel cuts through the IVC, separated by 0.7 arc- 
sec, and use these to make up an RA-Dec rectangle of size 
~ 5" x 20" within which cloudlet sizes could be obtained. 
Unfortunately, due to the presence of Sahara dust and non- 
ideal seeing conditions (1.5"-2.0"), integration times were 
somewhat longer than had been hoped for and in the end 
only two longslit positions were obtained, offset by a posi- 
tion angle of 85 degrees. Additionally, on the second night 
we obtained 1200-s of integration towards a blank sky region 
with a 0.85" slit and on the third night a total of 5x400- 
s, interleaved with the object exposures and using a slit of 
width 4.0", corresponding to ~9500-s of integration with a 
slit of 0.85" used for the IVC observations. 

Data reduction was performed using standard meth- 
ods within irafQ. The aims were firstly to obtain a spec- 
trum of the inner part of the globular cluster over the whole 
wavelength region, and more importantly, produce longslit 
spectra of the Call lines alone for the two longslit positions 
taken. Data reduction included debiasing, and use of the sky 
flatfields to get rid of the vignetting that was clearly visi- 
ble on many of the orders. Cosmic ray hits near the Call 
line were removed within FIGARO (Shortridge et al. 1999) 
using CLEAN. Because of the small pixel-to-pixel variation 
and pixel size (~ 0.2"), no flatfielding to remove the pixel- 
to-pixel variation was performed on the data. The images 
were rotated in order to make the cross-dispersion axis oc- 
cur along the image rows, after which sky subtraction was 
performed. As expected, this proved to be challenging. 

2.3.1 Sky subtraction 

Two methods were used to estimate the sky value and the re- 
sults compared. The first used the data from the third nights 
observation where the blank sky exposures were interleaved 
with the object exposures. For these data, we scaled the high 
signal-to-noise combined sky flat (taken at twilight with the 
0.85 arcsec slit) to produce an image 'scaled_sky' and then 
smoothed this in wavelength to a resolution equivalent to 
a 4 arcsecond slit. The scaling factor was chosen so that 
the resulting image was the same as the blank sky image 
taken with the 4 arcsecond slit for an equivalent of 5400-s 
of integration. The image 'scaled_sky' was then subtracted 
from each 5400-s of the object data. This method has the 
advantage that in the blank sky images there are clearly no 
contaminating counts from the object, but is susceptible to 
changes in the observing conditions and requires careful re- 
moval of the bias strip in both sky and object frames. Fig. ^ 
shows a cut across the dispersion axis for the sky exposures 
on night 3. As can be seen, there are few sky counts, al- 
though obviously in the outer parts of the slit there are few 
object counts too, so errors in the equivalent width increase 
rapidly in the extracted spectra for this region. For an equiv- 
alent 5400-s of integration time, at A ~ 3937A, there were 
~ 3 ADU sky counts and somewhat less at A=3933.66A, 

t iraf is distributed by the National Optical Astronomy Obser- 
vatories, U.S.A. 
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Figure 2. Cut across the dispersion axis at ~ 3937A of a blank 
sky region, equivalent to 5400-s of object observing time on night 
3. The 'dips' in the spectrum correspond to the interorder gaps. 

due to the 'sky' effectively being a solar (G-type) spectrum 
which has strong Call absorption lines. 

The method of sky subtraction that we had envisaged 
performing a priori involved removal of sky using parts of 
the longslit that appeared to be free of emission from the 
globular cluster. The problem with this method is illustrated 
by Fig. ^|(a) , which shows cuts across the dispersion axis for 
5400-s of data taken on the third night at each of the two 
slit positions. As is apparent from the Figures, the 'sky' level 
at the same position on the order at the two different slit 
positions is different by ~ 3-4 ADU per 5400-s of integra- 
tion. During night 3, at A ~ 3937A, position pi had ~ 8 
ADU as its 'sky' value, with the corresponding value for p2 
being ~ 12 ADU. These clearly are much larger than the 
values obtained by exposing on the night sky and combined 
with the fact that pi and p2 are different indicates that we 
are still obtaining counts from the cluster in the outer parts 
of the slit. This is confirmed by performing an extraction 
using such counts; in some parts of the spectrum the IVC 
absorption becomes negative which is unphysical. 

Because of the lack of object-free emission on the slit, 
we decided to remove the sky by subtracting the smoothed 
version of the sky images from our object data. This method 
relies upon the sky conditions being similar throughout the 
three nights of observing. This was checked by taking cuts 
across three places on the dispersion axis of each 5400-s 
worth of data and checking for variability. At A ~ 3937 A, 
the values varied by only 2 ADU per 5400-s between nights 
1 and 3 for the two slit positions. This variability, combined 
with the error in removing the bias, gives us an error in the 
final sky value of ~ 4 ADU at 3937A, compared with a peak 
continuum value at this point of ~ 30 ADU. 

2.3.2 Extraction and Analysis 

The sky-subtracted spectra were extracted within FIGARO by 
summing over each 10 columns, corresponding to the worst 
seeing of ~ 2.0 arcseconds. Wavelength calibration was then 
undertaken, taking into account the shift in the A scale along 
the length of the slit. Finally, normalisation and profile fit- 
ting were performed within DIPSO (Howarth et al. 1996). 
The final spectra have spatial resolution of ~ 2 arcseconds 
and velocity resolution of ~ 7 kms -1 . 
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Figure 3. (a) Cut across the dispersion axis at ~ 3937A for the 
two slit positions with 5400-s of integration, showing the difficulty 
in determining the sky. IOG refers to the interorder gap between 
the echclle orders (b) Cut along the dispersion axis before nor- 
malisation. 



Data were analysed using standard profile fitting meth- 
ods within DIPSO using the elf and IS programs. Before 
fitting, the the stellar line was removed by fitting the pro- 
file by eye to the lower-wavelength part of the spectrum 
up to the peak stellar absorption unaffected by LV and IV 
gas, and creating a mirror of this spectrum for the higher 
wavelength data. The derived stellar spectrum has a stark- 
like line profile, typical of stellar Call lines, and was sub- 
tracted from the whole spectrum to leave the interstellar 
components only. After stellar-line removal, the elf rou- 
tines were used to fit Gaussian profiles to input spectra and 
hence provide the equivalent width, peak absorption and 
full width half maximum velocity values for the low velocity 
gas and IVC components. These results were compared to 
the theoretical absorption profiles computed by the IS suite 
of programs within DIPSO. These theoretical profiles are de- 
rived from the observed b- values (corrected for instrumental 
broadening), atomic data for the Call lines taken from Mor- 
ton (1991) and initial guesses for the Call column density. 
Input parameters were varied until a good fit was produced 
to match the profiles obtained using elf. This method is 
only appropriate in the regime where the lines are not sat- 
urated and are resolved in velocity; both caveats apply for 
a number of positions in the current dataset, although the 
presence of narrow, unresolved components cannot be ruled 
out. If present, these would cause our estimated column den- 
sities to be too low. The final elf and is fits provide the 
Call number densities and values for the FWHM velocities 
at each of the positions along the slit. 



2.4 Wisconsin Ha Mapper observations 

Data were retrieved from the Wisconsin Ha Mapper 
(WHAM) survey in the range +06° < Dec. < +15° and 
2 Qh 2Q m < RA < 2i h 38 m (J2000). These data are of reso- 
lution 1 degree spatially and ~ 12 kms" 1 in velocity, with 
a velocity coverage at this sky position from ~ —120 to 
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Figure 4. WSRT Hi column density map of the M 15 IVC 

integrated between +50 and +90 kms -1 and at a resolu- 
tion of lll"x56" (~ NS by EW). Contour levels are at 
N H /=(2,4,6,8,10)xl0 19 cm" 2 . 

+90 kms -1 in the LSR. Data were reduced using standard 
methods, which included conversion into the LSR and re- 
moval of the geocoronal Ha line that appears at velocities 
in the range ~+30 to +40 kms" 1 (Haffner 1999; Haffner et 
al. 2001b). Due to the removal of this line, the noise in this 
part of the spectrum is slightly enhanced compared with the 
typical RMS noise value (measured for the current spectra) 
of ~1— 3 mR (kms" 1 )" 1 . These values are typical for the 
WHAM survey as a whole. Finally, we note if the IVC Ha 
flux originates from above the Galactic plane, then a cor- 
rection for Galactic extinction of 1.25 would be necessary, 
calculated from a E(B — V)=0.10 observed towards M 15 by 
Durrell & Harris (1993), combined with the extinction law 
of Cardelli, Clayton & Mathis (1989). We have not applied 
this correction to our data so the IVC Ha fluxes may be ~ 
25 per cent too low. 



2.5 IRAS ISSA data retrieval 

Both the 60 and 100 micron images of a field of size 
0.675° x 0.625° were extracted from the on-line versions of 
the IRAS Sky Survey Atlas (ISSA). These images are of res- 
olution ~ 5 arcmin and have most of their zodiacal emission 
removed although they have an arbitrary zero-point. 



3 RESULTS 

3.1 WSRT aperture synthesis Hi results 

Fig. ^ shows the H I column density map of the IVC for the 
low-resolution WSRT data alone, with Fig. ^| showing the 
corresponding map of the combined (WSRT plus Arecibo) 
image overlaid on the digitised Palomar Observatory Sky 
Atlas (POSS-I) red image regridded to a 4 arcsec pixel size. 
The major axis of this part of the IVC lies parallel with 
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Figure 5. WSRT plus Arccibo combined Hi column density map 
of the M 15 IVC integrated between +50 and +90 kms" 1 and 
at a resolution of lll"x56", overlaid on the red POSS-1 digital 
sky survey image towards M 15. Velocities are in the LSR with 
contour levels being at N HJ = (2,4,6,8, 10,12, 14) xlO 19 cm~ 2 . Lines 
of constant Galactic longitude and latitude are also plotted on the 
figure. 



the Galactic plane. The Hi channel maps of the combined 
WSRT plus Arecibo dataset in brightness temperature (Tb) 
are shown in Fig. [], where flux density per beam is related 
to Tb by S mJy/bcam =0.65fi as T B /A 2 m ; here fi as is the beam 
area in arcsec 2 and A 2 m is the observed wavelength in cm 
(Braun & Burton 2000). Immediately obvious from each of 
these figures is the fact that the H I is clumpy in nature, as is 
seen in other intermediate and high velocity clouds observed 
at a similar resolution (e.g. Wei/3 et al. 1999; Braun & Bur- 
ton 2000). Variations in the column density of a factor of ~ 4 
on scales of ~ 5 arcmin are observed, corresponding to scales 
of ~ 1.5 D pc, where D is the the IVC distance in kpc. The 
peak brightness temperature in the combined map (of reso- 
lution lll"x56"), is ~ 8 K, rising to ~ 14 K for the highest 
resolution WSRT map with beamsize 81"xl4". These val- 
ues compare with peak values of Tb for HVCs of ~ 25 K 
towards Complex A at 2 arcmin resolution (Schwarz, Sulli- 
van & Hulsbosch 1976), ~ 34 K observed towards Complex 
M at a resolution of 1 arcmin (Wakker & Schwarz 1991), 
and ~ 75 K observed to the compact high velocity cloud 
CHVC 125+41-207 (Braun & Burton 2000). Similarly the 
IVC observed by Wei/3 et al. (1999) has peak Tb exceeding 
7 K, at a resolution of 1 arcmin. 

The peak Hi column density in the combined image is 
~ 1.5x 10 20 cm -2 . If we assume that the cloud is spherically 
symmetric or filamentary (Stoppelenburg, Schwarz & van 
Woerden 1998), and using a cloud size full width half maxi- 
mum (c.f. Wakker & Schwarz 1991) of ~ 7 arcmin or 2 Dpc 
for the brightest point, we obtain a peak volume density of 
~ 24D -1 cm -3 . At the assumed distance, this value is a 
lower limit as it is likely that there will be more structure 



on smaller scales, as indicated by the spectra of Meyer & 
Lauroesch (1999). 

Fig. ^ shows the brightness temperature (Tb) profile at 
the position of peak temperature in the uniformly-weighted 
data, with Fig. |§| depicting the combined temperature pro- 
files at a number of positions in the cloud. The full width 
half maximum (FWHM) varies from 5 kms -1 at the posi- 
tion of the three peaks in brightness temperature (Fig. and 
Fig. ^(a)-(c)), to more than 12 kms -1 at other locations in 
the IVC (Fig. |(d)). Some parts of the cloud (cf Fig. |(b)) 
are well-fitted by a Gaussian, whereas others, such as Figs. 
^|(e), appear to be made up of two narrowish components. 
Although the signal-to-noise in this latter part of the cloud 
is low, both the WSRT and Arecibo datasets show evidence 
for a two-component velocity substructure, perhaps indica- 
tive of H I cloudlets or of overlapping clouds in the line of 
sight. 

Fig. displays the results of single-component Gaussian 
fitting to the IVC using XGAUSS and shows how the smaller 
values of FWHM (and hence implied kinetic temperature) 
tend to occur where the H I is brightest in the three cores. At 
these positions, in the absence of turbulence, the kinematic 
temperature is of the order of 500 K. The M 15 IVC does 
not show any regular velocity gradients across its field as 
has been seen in many other HVCs and IVCs (e.g. Wei/3 et 
al. 1999; Brims et al. 2000). 

We finally note that an interesting feature of the H I 
surface density map is that the IVC is approximately cen- 
tred upon the globular cluster, with the regions of lowest 
column density being located towards the globular cluster 
centre. It is tempting to suggest some mechanism whereby 
the IVC and globular cluster are associated, for example by 
capture of the IVC by M 15, with the gas in the centre being 
ionised by the cluster itself. However, the existence of other 
IVC components in the nearby field, combined with the fact 
that M 15 is located at a radial velocity of ~ -100 kms -1 
(Harris 1996), compared with the the IVC at ~ +70 kms -1 , 
makes it likely that the two objects are simply line-of-sight 
companions. 

3.2 Lovell telescope multibeam Hi results 

The Hi column density map of a 4.9x9.3 degree field, with 
spatial resolution 12 arcmin, centred to the South- West of 
M 15, integrated from +50 to +90 kms -1 , is shown in Fig. 

No new firm Hi detections were obtained that are not 
already given in the Leiden/Dwingeloo survey (Hartmann 
& Burton 1997), Kennedy et al. (1998) or Smoker et al. 
(2001a). However, the extent of an IVC component around 
RA=21 h 27 m , Dec=+09°10' (J2000; feature 'D') was deter- 
mined to be ~ 0.7 degrees in RA-Dec at full width half 
power in the Nhi map. An example Hi spectrum in this 
area is shown in Fig. [j]. At this position, a single-component 
Gaussian fit within DIPSO gives a velocity centroid of 68+0.5 
kms -1 , FWHM of 12.0+1.2 kms -1 and peak T B of ~ 1 K, 
similar properties to the main complex studied in this paper. 
The peak IVC H I column density towards this component 
('D') is ~ 3xl0 19 cm -2 . 

Also plotted on Fig. [jl] are 'A', 'AR' and 'HD'. 'A' corre- 
sponds to the main clump observed by Kennedy et al. (1998) 
which has FWHM ~ 12-15 kms -1 at 12 arcmin resolution. 
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Figure 6. WSRT plus Arecibo combined Hi brightness temperature channel maps towards the M 15 IVC. Velocities are in the LSR and 
the velocity resolution is 1.03 kms" 1 . Contour levels are at T s = (-l,l,2,4,8) Kelvin. The spatial resolution is lll"x54" (~ NS by EW). 



The peak IVC column H I density derived using the current 
observations towards component 'A' of 3.7xl0 19 cm -2 is 
close to the value of 3.9xl0 19 cm -2 observed by Kennedy 
et al. (1998) with the same telescope and gives us faith in 
our calibration. 'AR' refers to previous IVC detections us- 
ing the Arecibo telescope which have FWHM=12 kms -1 
and vlsr=+(>1 kms -1 at a resolution of 3 arcmin. 'HD' 
is the position of the halo star HD 203664 (of distance 3 
kpc) in which two strong interstellar IVC Call K absorp- 
tion components are seen with FWHM=2.8 and 3.2 kms -1 
and v~ +80 and +75 kms -1 respectively (Ryans, Sembach 



& Keenan 1996). As this is an absorption-line measurement 
towards a star, the resolution is sub-arcsecond. Finally, we 
note that there is a hint of emission at RA=21 h 23 m 04", 
Dec=+14°12'42" (feature 'E' on Fig. |Io|) although this is 
very close to the noise and may be a spurious detection. 



3.3 William Herschel Telescope longslit UES 
results 

Equivalent widths of some strong stellar lines from the inner 
3 arcseconds of the first slit position are shown in Table LI 
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Figure 7. WSRT uniformly- weighted H I brightness temperature 
profile of the M 15 IVC at the position of peak Tb- The beamsize 
is 81"xl4". The dotted line depicts the Gaussian fit to the spec- 
trum and has FWHM=5.2±0.2 kms -1 . The dashed line shows 
the (data— single Gaussian fit) residual. 



Table 1 . Equivalent widths of some strong lines towards the cen- 
tre of M 15. Where there are two values present, these refer to 
results taken on nights 1 and 2 respectively. Stellar types cor- 
responding to the EWs for that particular line are taken from 
Jaschek & Jaschek (1987) (JJ87). 



Species 

Can K 
Fel 
Cai 



■^rest 

3933.66 
4045.82 
4226 



EW (A) JJ87 Type 



1.6, 
0.16, 
0.25, 



1.5 

0.17 

0.28 



A2 
A3 
F2 



As the core of the cluster is unresolved from the ground, 
the spectrum obtained is composite. Fig. [n] shows the Can 
K spectra at this position. The equivalent widths of the LV 
and IV interstellar components at this position are ~ 0.3 
A and ~ 0.08 A respectively. The strength of the LV com- 
ponent multiplied by sin(b) of ~ 135 mAcompares with the 
canonical value of ~ 110 mA integrated EW of the Call K 
perpendicular to the Galactic plane (Bowen 1991). 

We now turn to the interstellar Ca n observations. Fig. 
|l2] shows the extracted spectra at three and five positions 
along the dispersion axis for slit positions 1 and 2 respec- 
tively, where the signal-to-noise in the continuum exceeds ~ 
10. The positions shown are separated by 2 arcseconds which 
corresponds to the worst seeing during the run, and is also 
the spatial resolution to which the data were smoothed when 
the extraction was performed. 

The prime aim of these longslit observations was to de- 
termine whether there is velocity substructure within the 
IVC. As can be seen from a few of the sightlines, there is 
tentative evidence for such structure, with two cloudlets be- 
ing present at LSR velocities of ~ +52 kms -1 and ~ +66 
kms" 1 (Fig. This corresponds to a separation of ~ 8 
pixels so is not an artefact caused by shifts in the echelle over 
periods of several hours. However, the signal-to-noise of the 
data is low, and follow-up observations are required to con- 
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Figure 8. WSRT plus Arecibo combined Hi brightness temper- 
ature profiles towards selected regions of the M 15 IVC. (a) First 
T s peak (RA=2l' t 29 m 44 s , Dec.=+12°09'42"). (b) Second T s 
peak (RA=21 h 29 m 36 s , Dec.=+12°06'22"). (c) Third T B peak 
(RA=21 h 30 m 13 s , Dec.=+12°10'42"). (d) High FWHM region 
(RA=2l' I 29 m 21 s , Dec.=+12°04'22"). (e) Two component pro- 
file (RA=2l' l 30 m 38 s , Dec.=+12°19'22"). The equinox is J2000. 



firm this finding. We note that although the observations 
of Meyer & Lauroesch (1999) may have also been expected 
to find such velocity substructure at their resolution of 14 
kms - , their observations were in the Nal D lines, which in 
some cases do not show velocity components, even when the 
Call lines do (Ryans et al. 1996). 

The current observations could also be used to estimate 
the IV C Cal l column densities along the slit as explained in 
2.3.2. Figure ^ displays the equivalent widths and 



section 

Call column densities for a number of the slit positions for 
the second slit orientation. Our Can K equivalent widths 
ranee from 0.06-0.09 A. These are slightly lower than those 
measured by Lehner et al. (1999) towards other parts of 
the IVC with low-resolution data, which are in the range 
0.05-0.20 A, with a median value of 0.10 A. Towards the 
nearby halo star HD 203664 (which has an IVC Hi col- 
umn density of less than ~ 10 18 cm -2 ), Ryans et al. (1996) 
measured an IVC Can equivalent width of 0.06+0.01 A. 
The Call/Hi ratio for the current data towards the M 15 
IVC varies from 2.1xl0~ 8 to 2.9xl0" 8 . Unfortunately, poor 
quality of our measurements caused by low signal-to-noise, 
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Figure 9. WSRT plus Arecibo Hi combined image. Greyscale; 
full width half maxima velocity widths obtained using single- 
component Gaussian fitting. Contours; peak values of bright- 
ness temperature of the same single-component fits with levels 
at (2,4,6,8) Kelvin. 
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and uncertainties in determining the sky and continuum lev- 
els, results in our errors being larger than the differences 
between these values. 



3.4 Wisconsin Ha Mapper (WHAM) Results 

3.4.1 The mam IVC towards M 15 

Of the 184 spectra extracted from the WHAM survey, a 
clearly-defined Ha component at intermediate velocities ex- 
ceeding +50 kms" 1 is only obvious at three positions (AO, 
A2, A5) with a weak detection towards (A7). We note that 
because M 15 itself is at an LSR velocity of ~ —100 kms -1 , 
it does not contaminate the observed spectrum. The spec- 
tra observed at these positions are shown in Fig. [l5](c). Fig. 
^(a) depicts the locations of these WHAM pointings rela- 
tive to the Lovell telescope H I surface density map, with in- 
tegration of velocities from ~ +50—90 kms" 1 in the WHAM 
Ha data also being superimposed. Fig. |l^(a) also shows the 
WHAM data integrated from ~ +50—90 kms -1 , over the 
whole field mapped by the multibeam observations. The gas 
with brightness ~ 0.1 R is quite extended about the main 
IVC condensation, with an additional weak signal detected 
towards the tentative Hi detection 'E' (see Fig. |Io| ). 

Results of three-component Gaussian fitting at the four 
positions are shown in Table ^ uncorrected for Galactic ex- 
tinction. These fits take into account the instrumental pro- 
file of the instrument, which is comprised of a 12 kms -1 
Gaussian with low-level wings superimposed. Although al- 
most equally-well fitted using just two components (for the 
LV and IV gas), we chose to fit using three components 
due to the slight asymmetry in the LV component, most 
clearly seen towards position (A7). Table also shows the 



Figure 10. Lovell telescope multibeam H I column density map of 
resolution 12 arcmin in the general direction of the M 15 globular 
cluster integrated from +50 to +90 kms -1 in the LSR. Contour 
levels are at (l,2,3)xl0 19 cm -2 , with greyscale leve ls fr om 0— 
4xl0 19 cm' 2 . For the meaning of the labels see Sect. 3.2. 
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Figure 11. Call K central 3" towards M 15. 
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Figure 12. Normalised spectra each 2 arcseconds along the slit. 
Spectra are offset in units of 1.0 for clarity, (a) Can K spectra at 
slit position 1. (b) Call K spectra at slit position 2. 
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Figure 13. Solid line: Call K LV and IV spectrum. Dotted line; 
ELF multicomponent fit. Dashed line; (data— model) residual plus 
1.0. A possible second intermediate velocity component is marked 
by IV?. 
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Figure 14. Ca II fit results along the second slit position. 



central velocity and velocity widths observed in H I towards 
the component (A5), obtained by smoothing the Lovell tele- 
scope multibeam data to a resolution of 1 degree. 

The spectrum with the strongest IV Ha emission (A5) 
is depicted in Fig. |l^. This is the nearest grid position to 
the M 15 IVC, whose centre lies approximately at Z=65.01°, 
b = —27.31°. The peak IV Ha brightness in this direction 
obtained using three-component Gaussian fitting is 0.035 
Rayleighs (km s" 1 )" 1 for the low velocity gas, and 0.033 
Rayleighs (km s" 1 )" 1 for the intermediate velocity gas, 
where 1 Rayleigh is 10 6 /47r photons cm -2 s" 1 sr" 1 . The inte- 
grated fluxes are 2.2±0.1 and 1.3±0.1 Rayleighs for the (to- 
tal) LV and IV gas respectively. The Gaussian fit gives cen- 
troids of -53.9±4.6, -4.4±0.7 and +64.3±0.4 kms" 1 , and 
FWHM velocity widths of 26.7±10.4, 47.8±2.6 and 31.5±1.4 
kms -1 , for the low and intermediate velocity gas respec- 
tively. The velocity centroid at this position (A5) agrees 
within the errors with the H I data smoothed to the WHAM 
spatial resolution; this contrasts with the results of Tufte, 
Reynolds & Haffner (1998) who tentatively found an offset 
in velocity between Ha and Hi velocities of ~ 10 km -1 to- 
wards various HVC complexes. For the position (A5), the 
velocity width of the IV Ha spectrum of ~ 32 kms" 1 is 
some 10 kms -1 greater than the Hi data at 1 degree res- 
olution. For a gas at ~ 10 4 K, some ~ 22 kms -1 of this 
is caused by thermal broadening, with the remaining width 
being due to non-thermal motions and beamsmearing of dif- 
ferent components. The difference between the H I and Ha 
widths may imply that the two phases are not mixed. How- 
ever, at least qualitatively, there is reasonable coincidence 
between the Ha and Hi peaks (Fig. [l5](a)), although the 
mapped area is small. 

The second positive detection (towards (A2) in Fig. |l5|) , 
occurs in a region where the local H I column density, at the 
Lovell telescope resolution of ~ 12 arcmin, is lower than 
lxlO 19 cm -2 . However, when smoothed to a WHAM reso- 
lution of 1 degree, the Hi column density at this point is ~ 
l.Ox 10 19 cm -2 . Finally, towards (A7) there is a weak detec- 
tion of ~ 0.25 R at ~ +59 kms" 1 . At the WHAM resolution, 
the H I column density at this point is ~ 0.4 x 10 19 cm -2 . We 
note that this position is close to the Arecibo-measured H I 
position 'AR' (Fig. hob , which has a similar velocity centroid 
of vls.r=+61 kms" . 
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Figure 15. (a) Lovell-telescope IV Hi column density map taken from Kennedy et al. (1998) in the vicinity of M 15 with WHAM- 
observed positions marked (AO) — (A7). Contour levels (full line) are at Nij/=(1,2,3) X 10 19 cm -2 . The dashed lines show the Ho contours 
from the WHAM Ha observations, integrated from ~ +50—90 kms -1 with levels at 0.1 X (1,2,3,4,5) R. (b) Ha brightness integrated 
from ~ +50—90 kms -1 for the whole field observed with the Lovell telescope in Hi. Contour levels are at 0.1 X (1,2,3,4,5) R. (c) WHAM 
Ha spectra at the positions (AO) — (A7). Velocities are in the LSR. Where Ha is observed (positions A0,A2,A5,A7), a three-component 
Gaussian fit is indicated by a dotted line with the (data— model) residual shown by the dashed line. 



3.4.2 WHAM pointings towards HD 203664 and position 
'D' 

The nearest WHAM Ha pointings in the vicinity of the halo 
star HD 203664 are some 0.38 and 0.66 degrees from the 
stellar position and are depicted in Fig. [l7| Although there is 
a slight excess of Ha with velocities exceeding ~ +50 km s -1 , 
this is very close to the noise and at lower velocities than 
the Call K absorption seen by Ryans et al. (1996) which 
had velocities of +75-80 kms -1 . Fig. |l| shows the WHAM 
pointings superimposed on the Lovell telescope multibeam 
Hi column density map for the IVC towards IVC position 
'D' (c.f. Fig. Juj). Again, although there is some Ha excess 
towards position (d3), this is very weak (brightness 0.07 R), 
and is also at a low- velocity of ~ +50 kms -1 ; this compares 
with the IVC Hi velocity of ~ +68 kms -1 at this point. 
Even though the WHAM pointing centre is just within the 
Nhj = 1x10 19 cm -2 contour, the relatively small size of the 



'D' component results in the H I column density at this point 
at 1 degree resolution being ~ 0.5 xlO 19 cm -2 . 



3.5 IRAS ISSA results 

In this section the IRAS 60 and 100 micron ISSA images are 
compared with previous Arecibo H I observations of Smoker 
et al. (2001a). These latter data are of resolution 3 arcmin 
and contain both low and intermediate velocity gas. Fig. |l9| 
shows the 100 micron image overlaid on the total (LV plus 
IV) Hi column density, with Fig. ^ displaying the corre- 
sponding 60 micron data and the IVC Hi column density 
alone. There appears to be a relatively good correlation be- 
tween the IRAS and H I results, both for the low and inter- 
mediate velocity gas; although with such a small field such 
agreement could have occurred by chance. The emission is 
not likely to be from the M 15 globular cluster itself, as glob- 
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Table 2. Wisconsin Ha Mapper results. The values for the velocities and fluxes are the results of three-component Gaussian fitting to 
the low and intermediate velocity gas. Ha brightnesses have not been corrected for Galactic absorption, hence assuming that the IVC 
lies above the Galactic plane, the tabulated values are likely to be low by ~ 25 per cent. Lovell telescope Hi multibeam results, smoothed 
to a WHAM-sized 1 degree beam, provide Hi velocity information towards position (A5). 



Position on Fig. [15j 

RA (J2000) 

Dec (J2000) 

I (degrees) 

b (degrees) 

Ha LV V c (kins" 1 ) 

Ha LV FWHM (kms" 1 ) 

Ha LV Flux (R) 

Ha LV Peak (mR (kms" 1 )- 1 ) 

Ha IV V c (kms" 1 ) 

Ha IV FWHM (kms" 1 ) 

Ha IV Flux (R) 

Ha IV Peak (mR (kms- 1 )- 1 ) 

Hi IV V c (kms- 1 ) 

Hi IV FWHM (kms" 1 ) 



(A7) 
21 h 26 m 53.3 s 
+11°28'58" 
63.88 
-27.16 

-32.0±16.6, +2.7±2.6 
44.6±22.4, 35.1±4.8 
0.27±0.18, 1.02±0.20 
4.8±1, 23±1 
+59.3±1.1 
18.8±3.4 
0.25±0.02 
11±1 



(A5) 
2l' I 29 m 24.9 s 
+12° 14' 30" 
64.98 
-27.16 

-53.9±4.6, -4.4±0.7 
26.7±10.4, 47.8±2.6 
0.14±0.06, 2.10±0.08 
4.2±1, 35±1 
+64.3±0.4 
31.5±1.4 
1.30±0.04 
33±1 
65.7±0.6 
22.0±2.0 



(A2) 
2l' i 32 m 00.4 s 
+13°00'02" 
66.09 
-27.16 

-43.8±3.7, -5.3±0.6 
18.8±9.4, 35.2±2.0 
0.11±0.04, 1.50±0.06 
4.7±01, 34±1 
+59.4±0.8 
38.1±2.4 
0.95±0.04 
20±1 



(AO) 
21 h 34 m 37.0 s 
+13°44'43" 
67.19 
-27.16 

-32.4±9.5, -0.7±1.5 
30.3±14.1, 29.5±2.9 
0.20±0.13, 1.26±0.14 
5.2±1, 34±1 
+71.4±0.5 
27.5±1.5 
0.91±0.03 
26±1 
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Figure 16. Solid line: WHAM Ha spectrum towards Z=64.98°, 
b=— 27.16° (A5). The dotted line shows the results of a three- 
component Gaussian fit. The dashed line is a (data-three compo- 
nent Gaussian fit) residual. 



ular clusters are undetected at the 60 and 100 micron IRAS 
sensitivity limits (Origlia, Ferraro & Pecci 1996). Plots com- 
paring the low, total and intermediate Hi column density 
with the IRAS flux density are displayed in Fig. |2lJ As- 
suming that the observed correlations are in fact real, the 
current data indicate that this IVC contains dust, as has also 
been observed in other such objects (Malawi 1989; Wei/3 et 
al. 1999). The fact that the IVC is tentatively detected at 
60 microns indicates that the dust may be warm, perhaps 
heated by collisions. 

The (tentative) IRAS detection contrasts with the situ- 
ation in HVCs, which at least at the IRAS sensitivity do not 
appear to emit at these wavelengths (Wakker & Boulanger 
1986). Such a difference indicates either differences in dust 



Figure 17. Nearest two WHAM Ha spectra towards HD 203664 
with their angular distances from the star indicated. 

content, dust temperature or environment (such as distance 
from the heating field of the Galactic plane). Hi absorption 
line spectroscopy towards a number of HVCs by Akeson & 
Blitz (1999) indicates that the fraction of cold gas is low in 
these objects. If so, then this would point to the Galactic- 
halo HVCs containing less dust than the IVCs and not dif- 
ferences in temperature. 

4 DISCUSSION 

4.1 The Ha properties of the M 15 IVC 

4-1.1 Estimates of the emission measure and Hll column 
density 

The detection of Ha towards the main M 15 IVC (position 
A5) with brightness of ~ 1.3 R compares to lower values of 
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Figure 18. Top panel: Lovell telescope multibcam IVC Hi sur- 
face density map of feature (D) (+50< V LSR <+80 kms" 1 ), with 
contour levels at N H/ =(1,2,3) X 10 19 cm" 2 with WHAM point- 
ings (dl) — (d5) marked. Bottom panel: WHAM Ha spectrum to- 
wards (d3). The dotted line is a three-component Gaussian fit 
with the dashed line showing the data-fit residual spectrum. 



Figure 19. Arecibo Hi and IRAS ISSA data. Greyscale: IRAS 
100 micron flux density. Contour: Total (LV plus IV) H I column 
density at N H/ = (56,60,64,68,72) X 10 19 cm" 2 . 
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between 0.06 to 0.20 R observed in a number of High Veloc- 
ity Clouds (Tufte et al. 1998). Similarly, many of the HVCs 
observed by Weiner, Vogel & Williams (1999) are very faint 
(<0.03 R), implying that, if they are photoionised by hard 
photons from the Galaxy, they are 20-60 kpc distant. The 
relatively high Ha emission from the M 15 IVC in turn ar- 
gues for a closer distance, assuming that shock ionisation is 
not the dominant factor. This is a big assumption, as detec- 
tions of Ha in the Magellanic Stream (distance 55 kpc) of 
are 0.20-0.37 Rayleigh (Bland-Hawthorn & Maloney 1999) 
are similar to the range of ~ 0.1—0.5 Rayleigh detected 
towards the IVC Complex K (Haffner, Reynolds & Tufte 
2001a) which has a distance bracket of 0.3—7.7 kpc. 

Assuming that extinction is negligible and that the 
cloud is optically thin, the observed Ha brightness of the 
IVC can be used to estimate both the column density of 



Figure 20. Arecibo Hi and IRAS ISSA data. Greyscale: IRAS 
60 micron flux density. Contour: IVC H I column density at 
Nffj=(2,4,6,8)xl0 19 cm~ 2 . Note how the peak Hi contour cor- 
responds with the IRAS peak. 

ionised hydrogen and the electron density. For an optically 
thin gas, the Ha surface brightness in Rayleigh, is given 
by; 

I a = 0.36 x EM x %/T x (1 — 0.37 x ln(T)), (1) 
where T is the gas temperature in units of 10 4 K and EM is 
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N m IVC (10 19 cm" 2 ) 



Figure 21. Arecibo Hi and IRAS ISSA data. The IRAS fluxes 
have an arbitrary zero-point, (a) IRAS 100 micron flux density 
plotted against low velocity gas Hi column density, (b) IRAS 
100 micron flux density plotted against total H I column density, 
(c) IRAS 60 micron flux density plotted against IVC H I column 
density alone. 

the emission measure ( J n 2 d£) in units of pc cm -6 (Simon- 
etti, Topasna & Dennison 1996). Here, n e is the electron 
density and integration over dl provides the thickness (L) of 
the ionised region. Assuming the temperature of the Hn is 
8,000 K, typical of the warm interstellar medium (Reynolds 
1985), the resulting emission measure towards position (A5) 
is 3.7 pc cm -6 . 

Assuming a distance to the main IVC condensation 
of 1 kpc, the cloud size (d) at a column density limit of 
Nhi =10 19 cm -2 is ~ 0.6° xO. 8°, corresponding to linear 
dimensions of ~ 10x14 pc. If the cloud is approximately 
spherically symmetric, this results in an estimated electron 
density (assuming a size of 12 pc) of yjFM/d ~ 0.6/ _1/2 
cm -3 and a column density of ionised hydrogen (=L x n e ) of 
~ 2xl0 19 / 1/2 cm' 2 , where / is the filling fraction of Hn in 
the IVC. Smoothed to a resolution of 1 degree (correspond- 
ing to the size of the WHAM beam), the Hi column density 
at this point derived from the data of Kennedy et al. (1998) 
is ~ 1.4x 10 19 cm -2 , hence the fractional ratio of H II to H I 
is of order lAxf 1 ' 2 , indicating that there is a substantial 
amount of ionised gas at intermediate velocities present in 
this sightline. 

For positions (A2) and (A7), and using the same cloud 
size as above, results in estimated Hn column densities of 
1.5xl0 19 cm -2 and 0.8xl0 19 cm -2 , corresponding to frac- 



tional H II to Hi values of order 1.5x/ 1//2 and 2.0x/ 1//2 re- 
spectively. These fractional Hn to Hi ratios towards the 
current IVC are somewhat higher than derived by Tufte et 
al. (1998) for a sample of HVCs which were calculated to be 
< 0.06 x/ 1/2 . However, recent work by Bluhm et al. (2001), 
in sightlines towards the Large Magellanic Cloud, used the 
relative underabundance of neutral oxygen to infer an ionisa- 
tion level in both an IVC and HVC of ~ 90 per cent. It would 
be useful to observe this cloud using the same methods as 
described in the current paper and compare the results. 

Towards (d3) there is no detection in Ha with velocities 
corresponding to the Hi value of ~ +68 kms" 1 . Assuming 
a cloud size of 8pc and upper limit to the Ha brightness of 
0.2 R, gives an estimated upper limit to the IVC H II column 
density and fractional Hn to Hi values of ~ 4xl0 18 cm -2 
and 0.7x/ 1/2 . Finally, we consider the halo star HD 203664, 
towards which the limiting column density of neutral hydro- 
gen is ~ 10 18 cm" 2 . If we take the upper limit to the Ha 
brightness of 0.2 R, and combine this with a cloud size L 
(in pc, where the Arecibo beam is ~ 1 pc FWHM at 1 kpc 
distance), a fractional Hn to Hi ratio of ~ 1.4xL x / 1//2 can 
be set by the current observations. 



4-1-2 Estimates of the ionising radiation field and 
electron density 

Following Tufte et al. (1998), if photoionisation is the dom- 
inant cause of Ha emission, then the incident Lyman con- 
tinuum flux Flc can be estimated thus, assuming case B 
recombination; 

Flc = 2.1 x 10 5 -^_ photons cm' 2 s" 1 , (2) 
0.1R 

-where I a is the Ha intensity in Rayleigh. For the M 15 IVC, 
equation |2j implies an incident flux (Flc) of 2.7xl0 6 pho- 
tons cm - s _1 . Hence if photoionisation is the main cause 
of Ha production, the derived Lyman-alpha continuum flux 
towards the main M 15 IVC condensation is a factor 6—22 
times higher than the implied incident flux towards the A, C 
and M HVCs observed by Tufte et al. (1998), and more than 
twice that observed towards the Complex K IVC by Haffner 
et al. (2001a). We recall that Complex A lies between 4—10 
kpc, with Complex C being some 5—25 kpc distant and the 
M 15 IVC being closer than 3 kpc. In the future, comparison 
of derived Lyman-alpha continuum fluxes for a larger sample 
of IVC and HVC sightlines with known distances may pro- 
vide information on the relative contributions of the Galactic 
and extragalactic ionising field. 

An alternative possibility is that the Ha is produced 
by shock ionisation, caused by interaction of the IVC with 
LV gas. This is a real possibility given the orientation of 
the IVC and the fact that its z-distance of less than 1 kpc 
puts it in the lower Galactic halo. Towards the nearby halo 
star HD 203664 in which IV absorption is seen, Sembach 
(1995) postulated that the dust grains in the IVC have 
been processed by such shocks which also currently produce 
the highly ionised species. For an ambient temperature of 
<3xl0 5 K, a cloud of velocity 50 kms - will be supersonic 
and hence shocks may be formed given the right conditions. 
Using the models of Raymond (1979), which are applica- 
ble for shocks with speed 50<V S <140 kms -1 , the face-on 
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Ha surface brightness (I ape r P ) can be related to the number 
density of the pre-shocked gas (no), thus; 

(I ap er P ) ~ 6.5 x n (T4/100) 17 R (3) 

Given that I aperp =1.3 R, and assuming a shock speed 
of 50 kms -1 , leads to an upper limit to no of 0.6 cm -3 . The 
value is an upper limit (for this shock speed velocity) as 
a non-perpendicular sightline will increase the observed I a 
(Tufte et al. 1998). Of course, given the fact that the trans- 
verse component of the velocity of the M 15 IVC is unknown, 
this value is very uncertain. Discriminating between shock 
and photoionisation is difficult, although further progress 
may be possible via measurements of appropriate emission 
line ratios (Tufte et al. 1998). 



4.2 Can number density towards the M 15 IVC 
and HD 203664 

Before discussing the Call K results, we note that calcium 
is depleted onto dust and is not the dominant ionisation 
species, hence the absolute metallicity of the M 15 IVC is 
uncertain and awaits high resolution UV observations. As 
emphasized by Ryans et al. (1997), differences in resolu- 
tion between the optical and radio data, combined with 
the Call K results only placing limits on the ion abun- 
dances, makes it important not to over-interpret the ob- 
served N(Caii)/N(Hi) ratios. 

With that caveat, and assuming that the current ob- 
servations do not miss any narrow-velocity components, the 
average IVC Ca n number density and ratio of the IVC com- 
pared to Hi towards the centre of M15 are logio(N(Can 
cm" 2 ) = 12.0 and N(Caii)/N(Hi) = 2.5xl0" 8 . The Hi col- 
umn density towards the M 15 centre is 4xl0 19 cm" 2 and 
was obtained using the combined WSRT plus Arecibo map 
of resolution lll"x56". For the nearby sightline towards the 
halo star HD 203664, we use the results of Ryans et al. (1996) 
for the total IVC Call K column density of ~lxl0 12 cm" 2 , 
combined with the upper limit to the Hi column density 
at a resolution of 3 arcmin towards HD 203664 of ~lxl0 18 
cm -2 (Smoker et al. 2001a), giving a much higher value of 
N(Caii)/N(Hi) > l.OxlO -6 . 

The current results compare with literature values of 
N(Caii)/N(Hi) ~ 2xl0 -9 (Hobbs 1974, 1976) for low ve- 
locity diffuse clouds and of N(Caii)/N(Hi) in the range 3- 
300 xlO -9 cm" 2 for the high latitude clouds studied by Al- 
bert et al. (1993). The N(Call)/N(Hl) ratios in IVCs and 
HVCs are thought to be higher than in low velocity gas 
due to the former having less dust onto which calcium is de- 
pleted. Thus the observed N(Call)/N(Hl) ratio of 2.5xl0" 8 
(or ~ 0.01 of the total solar calcium abundance) in the line- 
of-sight towards M 15 is typical of other high latitude clouds 
and also of other HVCs and IVCs (e.g. Wakker et al. 1996a; 
Ryans et al. 1997). 

The lower limit of N(Caii)/N(Hi) > 10" 6 towards HD 
203664 is, however, on the high side for IVCs/HVCs, be- 
ing ~ 0.5 of the total solar calcium abundance. There are 
several possible reasons for this. Firstly, the (currently un- 
detected) Hi towards HD 203664 could be in a clump of 
gas smaller than the Arecibo beamsize of 3 arcmin; if this 
were the case then the H I column density limit used would 
be too low and the derived N(Call)/N(Hl) ratio too high. 



Additionally, HST UV observations towards HD 203664 in- 
dicate that the Hi towards this object is at least partially 
ionised (Sembach, private communication), either by shock 
ionisation or photoionisation. If photoionisation, aside from 
the normal ionising source being Galactic OB-type stars or 
the extragalactic UV field, HD 203664 itself (spectral type 
B0. 5) could be a possible ionising source. The fact that its 
LSR velocity is +110 kms -1 (Little et al. 1994) compared 
with the IVC at +70 kms" 1 is inconclusive in determining 
the relative distance of the line of sight IVC towards HD 
203664 with the star itself. Finally, there remains the pos- 
sibility that the hydrogen is in molecular form. However, 
given the low H I column density towards the HD 203664 
sightline, this appears unlikely. 

Alternatively, it could be that the derived value of 
N(Caii)/N(Hi) > 10" 6 towards HD 203664 is correct. This 
would tally with the IUE results of Sembach (1995), which 
found that the majority of the elements in the IV gas, when 
referenced to sulphur, were within a factor 5 of their solar 
values and strongly point to a Galactic origin for this part of 
the IVC. The fact that our derived value for N(Caii)/N(Hi) 
towards the M 15 IVC of 0.01 solar is much lower than to- 
wards HD 203664 is likely to be caused by different ionisa- 
tion fractions and dust contents, and/or differing formation 
mechanisms. Clearly the latter is speculative and requires 
follow-up high resolution UV observations towards the M 15 
IVC to determine the abundances of elements such as sul- 
phur and zinc that are not depleted onto dust. 

4.3 Velocity widths and temperatures towards the 
M 15 IVC and HD 203664 

Towards the main M 15 IVC condensation, (feature 'A' on 
Fig. |l(j) values of the Hi FWHM velocity width at reso- 
lutions of ~2xl arcmin range from 5—14 kms - , corre- 
sponding to maximum kinetic temperatures in the range 
~500-4000 K. Mid-way between the M 15 IVC and HD 
203664, the FWHM equals 12 kms" 1 at a resolution of 3 
arcmin, which corresponds to T*, ~ 3000 K (feature 'AR' 
on Fig. |l^). The current observations have additionally ob- 
served feature 'D' at a resolution 12 arcmin, with FWHM 
velocity width also of 12 kms" 1 , indicating gas of similar 
temperature. We note that each of these temperatures will 
be upper limits due to beamsmearing and turbulent veloc- 
ity components. Finally, towards HD 203664, Ryans et al. 
(1996) found cloudlets with FWHM of 2.8 and 3.2 kms" 1 in 
Call K, corresponding to upper limits for the kinetic tem- 
peratures of ~ 8000—10,000 K. Towards the same star, Sem- 
bach (1995) used the relative abundances of low-ionisation 
species to derive a temperature for the HD 203664 IVC of 
5300-6100 K. 

The higher IVC gas temperatures towards HD 203664 
than towards the M 15 IVC, feature 'D' or the intermediate 
position 'AR' could be interpreted as being caused by the 
former cloud being nearer to the ionising field of the Galac- 
tic plane than the other two IVCs (Lehner 2000). It seems 
more likely that the lower temperatures seen towards parts 
of 'A' and 'AR' are simply caused by shielding of parts of 
these gas clouds; shielding that is not possible towards the 
HD 203664 IVC because of the lower gas density there. A 
two-phase core-envelope structure for halo HVCs has often 
been proposed within the Galactic corona for 1< z <5 kpc 
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(e.g. Wolfire et al. 1995), where the two components of 100 
and 10,000 K are entrained by pressure from the hot Galac- 
tic corona. We note that the high-end Hi temperatures ob- 
served towards the M 15 IVC indeed occur towards its outer 
parts where the Hi column density is low and the FWHM 
velocity widths are uncertain. 



4.4 Comparison of H I properties with previous 
Galactic halo IVCs and HVCs 

In this section we compare the high-resolution Hi proper- 
ties of the M 15 IVC with other IVCs and HVCs known 
to lie within the Galactic halo and observed at compara- 
ble resolution. We note that it is likely that there are many 
different types of HVC, with recent work, for example, in- 
dicating that a number of the compact HVCs lie at dis- 
tances of several tens of kpc (Braun & Burton 2000). A 
literature search found the following objects with known dis- 
tances and observed in Hi at high resolution: the M 13 IVC 
(Shaw et al. 1996), the 4-Lac HVC 100-7+100 (Stoppelen- 
burg et al. 1998), the Low-Latitude Intermediate Velocity 
Arch (Wakker et al. 1996b), IVC 135+54-45 (Wei/3 et al. 
1999), HVC 132+23-211 (within Complex A; Schwarz et 
al. 1976) and the M92 HVC (within Complex C; Smoker et 
al. 2001b). 

Table ^| compares the properties of these IVC and HVC 
H I gas clouds known to exist in the Galactic halo. Inspect- 
ing the table, there are no clear differences in the H I prop- 
erties (column density, peak temperature) of the two types 
of objects, which show a large scatter both within IVCs and 
HVCs. Similarly, the velocity widths of all of the objects, 
barring the peculiar HVC100— 7 and the Mirabel cloud, all 
show minimum values of FWHM ~ 5 kms" 1 and indicat- 
ing that gas of temperature less than ~ 500 K is common in 
these objects. This is in contrast to the situation observed at 
lower resolution for some northern clouds, where IVCs tend 
to have smaller velocity widths than their HVC counterparts 
(Davies, Buhl & Jafolla 1976). 

If some IVCs and HVCs are formed from infalling gas, 
sweeping up high-z Hi as they fall towards the plane, or 
if they are objects within the plane, or if they are formed 
within a Galactic fountain, it seems plausible that they could 
be preferentially aligned with the Galactic plane. Of course, 
it must be taken into consideration that at high resolution, 
only small parts of cloud complexes are studied, and by 
chance, some of these components will be aligned with the 
plane. In any case, of the four objects in the sample with a 
well-defined axis and at high Galactic latitude, three have 
their major axis near-parallel with the plane. Although there 
exist a number of such objects observed at lower resolution 
with this orientation, to our knowledge no systematic survey 
has been performed determining the orientation parameters 
of HVCs. If performed, this could act as a further discrim- 
inator between HVCs known to exist in the Galactic halo, 
and the sample of HVCs postulated to lie at extragalactic 
distances. 

Summarising, at present there are too few high- 
resolution observations of Galactic halo IVCs and HVCs 
to determine differences in Hi properties and any relation- 
ship between the two types of object. However, as previously 
noted, the IRAS and Ha properties do appear to differ, al- 



though the number of objects studied in all three wavebands 
remains small. 



5 SUMMARY AND CONCLUSIONS 

The current Hi WSRT synthesis observations have shown 
that on scales down to ~ 1 arcmin, the M 15 IVC shows sub- 
structure, with variations in the column density of a factor 
of ~ 4 on scales of ~ 5 arcmin being observed, correspond- 
ing to scales of ~ 1.5 D pc, where D is the the IVC distance 
in kpc. Of course, this is not an unexpected finding, but 
once again demonstrates that great care must be taken in 
interpreting quantities such as cloud metallicities which are 
derived from a combination of low-resolution radio plus op- 
tical data. The Lovell telescope H I observations towards this 
cloud demonstrated how relatively large areas of sky can be 
mapped with the multibeam system in a short period of time 
in the search for IVCs and HVCs. These data showed that 
the M15 IVC has components spread out over several square 
degrees, with component 'D' being mapped for the first time 
at medium resolution (12 arcmin) and having a similar col- 
umn density to the IV gas centred upon M15 itself. Both 
the Hi emission-line and Call absorption- line data showed 
tentative evidence for velocity substructure, perhaps indica- 
tive of cloudlets. The Call/Hi value of ~ 2.5xl0 -8 towards 
the main M 15 condensation is similar to that previously ob- 
served in other IVCs and HVCs. Towards HD 203664, the 
observed lower limit of 10~ 6 is somewhat higher, although 
this may be caused by factors such as the Hi beam being 
unfilled or partial ionisation of the gas on this sightline. The 
Hi properties of the M 15 IVC are indistinguishable from 
HVCs, although with the lack of distance information to- 
wards most HVCs, comparisons are difficult. 

The tentative detection of infrared emission from the 
M15 IVC, as in other IVCs, does distinguish it from HVCs, 
and either points to the M15 IVC containing more dust, 
and/or being closer to the heating field of the Galactic plane 
than HVCs, which as a class of objects are not detected 
in the IRAS wavebands. Similarly, the relatively strong Ha 
emission (exceeding 1 Rayleigh) towards parts of the M15 
IVC, if caused by photoionization, may place it closer to 
the Galaxy than HVCs. Again, however, this finding is un- 
certain due to the problem in distinguishing photoionisation 
from shock ionisation, uncertainties in dust content, and dif- 
ferences in Hi volume densities in different objects studied 
thus far. 

Future work towards this cloud should include higher- 
signal-to-noise observations in the Call line in order to de- 
termine if the cloud velocity substructure tentatively found 
in the current observations is in fact real, and whether the 
Call/Hi ratio determined by the current observations is 
lower than towards the HD 203664 sightline. This should 
be combined with 12 CO(1-0) sub-mm observations in order 
to determine if molecular material exists towards the peaks 
in H I column density and out of which stars may form. The 
determination of the falloff in Hi column density, of the 
cloud to low column density limits would also indicate the 
ionisation properties of the object and whether or not there 
is any interaction between the M 15 IVC gas and low veloc- 
ity material. Finally, UV observations towards M 15 globular 
cluster stars, although difficult, would provide important in- 
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Table 3. Previous high-resolution Hi observations of IVCs and HVCs known to exist in the halo, compared with the current results. 
T p and Nm are the peak temperatures in K and peak Hi column densities in units of 10 19 cm -2 respectively. AVI and AV2 are the 
broad and narrow velocity widths in kms . The field "Parallel to the Gal. plane' only refers to the object observed at high-resolution 
and not to the whole complex in the case of the large HVCs. References; RKSD97, Ryans et al. (1997); SSH76, Schwarz ct al. (1976); 
SBK96, Shaw ct al. (1996); SSW98, Stoppelenburg ct al. (1998); SKD01, Smoker et al. (2001b); WHS96, Wakker et al. (1996b); WPS99, 
van Woerdcn ct al. (1999); WWHM99 Wei/3 ct al. (1999). GSR velocities are calculated via Y G SR=^LSR + 250xsin(«)xcos(6) kms" 1 . 
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formation on the absolute metallicity of the gas towards this 
object for comparison with the HD 203664 sightline. 
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